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SWING-TWIST DECOMPOSITION IN CLIFFORD ALGEBRA 


PRZEMYSLAW DOBROWOLSKI 


Abstract. The swing-twist decomposition is a standard routine in motion 
planning for humanoid limbs. In this paper the decomposition formulas are 
derived and discussed in terms of Clifford algebra. With the decomposition one 
can express an arbitrary spinor as a product of a twist-free spinor and a swing- 
free spinor (or vice-versa) in 3-dimensional Euclidean space. It is shown that in 
the derived decomposition formula the twist factor is a generalized projection 
of a spinor onto a vector in Clifford algebra. As a practical application of 
the introduced theory an optimized decomposition algorithm is proposed. It 
favourably compares to existing swing-twist decomposition implementations. 


1. Introduction 

Swing-twist decomposition of rotations is commonly used in context of humanoid 
motion planning. Consider movement of an arm reaching some predefined position. 
In order to displace it properly, a controller calculates twist factor of the corre¬ 
sponding rotation. Having this factor computed, a controller is then able to apply 
some corrections to the motion so that the arm is not unnaturally twisted during 
the motion. Swing-twist decomposition is an inherent part of a correction al¬ 
gorithm. It allows one to decompose an arbitrary rotation into a swing part (tilt 
of a given axis) and a twist part (rotation around a given axis). In this paper the 
decomposition is derived and discussed in terms of Clifford algebra. 

Swing-twist decomposition has already been considered for a few decades. Differ¬ 
ent authors have obtained equivalent formulas, in particular for quaternion algebra. 
Current literature tends to neglect a deeper consideration on the spin-twist decom¬ 
position. There is wide range of publications on humanoid motion planning. Most 
of these works relate to swing-twist decomposition in some way. In the paper, the 
most revelant approaches to swing-twist decomposition are quoted and compared. 
Starting from the recent PhD dissertation by Huyghe (JT], 2011), the swing-twist 
decomposition is introduced in quaterion algebra by a proposed therein projection 
operator. Unfortunately, the origin of the projection operator is not enough ex¬ 
plained. In two previous papers, Baerlocher ([2], 2001) and Baerlocher, Boulic (0], 
2000) investigate joint boundaries for ball-and-socket joints using swing-twist de¬ 
composition. Grassia ([4j, 1998) discusses features of a swing-twist decomposition 
in terms of an exponential map. However, the author does not present any re¬ 
lated formula for the decomposition. A classic work by Korein ([5], 1984) contains 
most of the initial results on body positioning and joint motion. Among others 
the work uses swing-twist decomposition of rotations to constrain movement of an 
elbow. 

None of these are generalized to Clifford or geometric algebras which has deeper 
consequences than previous results. In this paper the swing-twist decomposition 
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is derived as an inverse of a formula expressing the set of rotations which move a 
given initial vector to cover a given terminal one. 

2. Existing and related solutions 

Swing-twist decomposition splits a given rotation in two parts: a swing part and 
a twist part. A schematic view of a limb rotating with a ball joint is presented in 
figure 





By zero twist reference vector one refers to a base vector with respect to 
which the swing-twist decomposition is performed. Usually, this is the bone of a 
rotating limb. 

There are several existing approaches. These solutions differ in terms of perfor¬ 
mance, complexity and result exactness. 

2.1. Direct method. For quaternions one can make the following argument. Let 
q £ Iff be a quaternion and v be a zero-twist reference vector. In case of twist-after¬ 
swing type of decomposition (see section [6]) the initial vector becomes the given v 
vector and the terminal vector is w = qvq*. Schematic view is presented in figure 

El 



Figure 2. A direct decomposition of a quaternion 

Swing quaternion can be calculated with an axis angle representation of the 
quaternion. Here, the axis is a normalized vector perpendicular to v and w and the 
angle is equal to the angle between v and w so: 

v x w 
n ||u x w\\ 
cos (a) = v ■ w 

q s = cos(a/2) + sin(a/2)(n x i + n y j + n z k) 
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Twist quaternion is calculated by the following inversion formula: 

q = qtq s =>qt = qq 7 1 

qt = q(cos(a/2) - sin(a/2)(n x i + n y j + n-k)) 

In this formulation both trigonometric and inverse trigonometric functions are 
used. It is a practical disadvantage since it requires computationally expensive 
functions which can also face some accuracy issues. 

2.2. Huyghe’s method. Huyghe uses swing-after-twist decomposition type. First 
a simplified formula is derived which is a decomposition of an arbitrary quaternion 
q into a product 

q = q s qt 

with respect to a constant Z-axis aligned zero-twist reference vector (called there 
” a twist axis”). With the following coordinates 

q = w + x'l + yi + zk 
q t =w t + z t k 
q s = w s + a; s i + y s j 

the author shows using some algebraic transformations that the decomposition is 

±10 

w t = 7 , 

V'iU + z 1 

±z 

Zt = , „ = - 

Vto 2 + z 2 

w s = w t w + z t z 
x a = w t x - z t y 
Vs = w t y + z t x 

It must be noted that some of the results were not rigorously stated. It particular, 
scenarios when some coefficients are equal to zero are not discussed separately which 
leads to possible division by zero. 

In the second part of the thesis, Huyghe derives a generalized formula for swing- 
twist decomposition. Assume that axis-angle representation of quaternions is 

q s = [w s ,u s ] = [cos(^),u s sin(|)] 

T T 

qt = [W t ,v t ] = [cos( —), tit sin(—)] 

Multiplying both quaternions one writes 

(1) q = q s qt = [w s w t - v t ■ v s , w s v t + w t v s + v s x v t \ 

Huyghe notes that the axes v s and Vt are perpendicular so formula (|T|) simplifies to 

(2) q = [w s w t , w s v t + w t v s + v s x v t ] 

Next, the author introduces a new quaternion q p (no origin is provided) which is 
”a projected version of the initial quaternion q onto the twist axis” an) 

q p = [w, (v ■ u t )u t \ 

= [w s w t , (w s v t ■ Ut + w t v s ■ u t + (v s X Vt) ■ u t )u t \ 

= [w s w t , wi s ||u t || 2 sin(^)u t ] 
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= [w s Wt,W s Vt] 

Normalization of quaternion q p gives 

% _ [w s w t ,w 8 v t ] _ w s [w t ,v t ) _ 

IMI V w s w t + H^IMI 2 WsV w t + IHI 2 

which in result turns out to be twist quaternion. Note that in the above equation 
Huyghe does not consider w s = 0 which is a drawback. Remaining swing quaternion 
is calculated from 

Qa = qq* t 

which completes the decomposition. 

3. Preliminaries 

The presented results relate to C l 3 := Cf 3i0 (K) - Clifford algebra of 3-dimensional 
real space. Herein, a spinor is given by s = a + be i 2 + ce 23 + de 3 i G Spin(3) and a 
vector is given by v = .xeo + ye 1 + ze 2 . A rotation of a vector is given by Clifford 
multiplication: 

v' = svs- 1 

for an arbitrary vector v and a spinor s. 

In this paper UN 1 denotes the inverse of a spinor. In case of Spin(3) it is equiv¬ 
alent to conjugation: 

s -1 = s = a — be i2 - ce 23 - cie 3 i 

A spinor can be written as a sum of its scalar and bivector parts: 

s = Ho + H2 

where [s]o G Ao anc ^ H 2 G /\ 2 K 3 . In several places the Hodge star * operator 
is used. In an orthonormal basis it is defined as 

*(ei A e 2 A • • • A e k ) = e fc+ i A e fc+2 A • • • A e n 

but in this paper it is used only in relation to the bivector part of a spinor. In this 
case the formula can be simplified to: 

*[s] 2 = —e i23 [s] 2 G f\R 3 
1 

A convenient notation is used for normalized vectors. For a given non-zero vector 

v A 0, the normalization function is defined as: 


where ||u|| = y/v ■ v = y/vv is the length of vector v. Note that it is impossible 
to define a normalized zero vector. Let s be a spinor and v be a vector in Ct 3 . 
The rotation svs -1 of the vector v by the spinor s gives a rotated vector v' = 
a/eo + y'e 1 + 2 /e 2 which is equal to: 

{ x' = (a 2 — b 2 + c 2 — d 2 )x + 2 y(ab + cd) + 2 z(bc — ad) 

y' = ( a 2 — b 2 — c 2 + d 2 )y + 2 x(cd — ab) + 2 z(bd + ac ) 

z’ = (a 2 +b 2 — c 2 — d 2 )z + 2x{bc + ad) + 2 y{bd — ac) 

This paper if organized as follows: first we obtain the complete set of spinors 
which do not rotate a given non-zero vector. Next we derive a spinor which rotates 
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a given initial vector to cover a given terminal vector. By combining these two 
formulas we derive a formula which represents an arbitrary rotation by initial and 
terminal vectors (swing factor) together with axis rotation (twist factor). Finally 
swing-twist representation is inversed for an arbitrary rotation and a unique swing- 
twist decomposition of a spinor is obtained. In conclusion some applications of the 
proposed decomposition are presented. 

4. The set of spinors which do not rotate a given vector 

In this section we assume that v is a non-zero vector. The set of all spinors 
s £ Spin(3) which do not rotate the given vector v will be called an invariant set 
of spinors for a given vector. The following theorem will be proved: 

Proposition 1 (The set of spinors which do not rotate a given vector). Let v be 
a non-zero vector in G(.^. The complete set of spinors s such that svs -1 = v is 
determined by: 

s = exp(ei 23 on(u)) 

for all a £ [0; 27r). For a given element a £ C ^3 exponential of a is defined as: 

exp(a) := £ £ 

k =0 

Proof. The proof will use the coordinate expansion of spinor rotation formula and 
the normalization identity giving the following set of equations: 

' (a 2 — b 2 + c 2 — d 2 — l)x + 2 y(ab + cd) + 2 z(bc — ad) = 0 

(a 2 — b 2 — c 2 + d 2 — 1 )y + 2 x(cd — ab) + 2 z(bd + ac) = 0 

(4) 

(a 2 + b 2 — <? — d 2 — 1 )z + 2x(bc + ad) + 2 y(bd — ac) = 0 
„ a 2 + b 2 + c 2 + d 2 - 1 = 0 

The key idea is to extract simple relations from [4] yet avoiding high order equations. 
Denote Qi as the left side of zth equation of [I] A valid solution ( a,b,c,d ) must 
satisfy the following equation: 

(5) xQi + yQ2 + zQs — ( x 2 + y 2 + z^)Q 4 = 0 
formula ([5]) can be expanded: 

(a 2 — b 2 + c 2 — d 2 — l)x 2 + 2 xy(ab + cd) + 2 xz(bc — ad) 

+ (a 2 — b 2 — c 2 + d 2 — 1 )y 2 + 2xy(cd — ab) + 2 yz(bd + ac) 

+ (a 2 + b 2 — c 2 — d 2 — 1 )z 2 + 2 xz(bc + ad) + 2yz{bd — ac) 

- {x 2 + y 2 + z 2 )(a 2 + b 2 + c 2 + d 2 - 1 ) = 0 

(a 2 — b 2 + c 2 — d 2 — 1 — a 2 — b 2 — c 2 — d 2 + l)x 2 + 2 xycd + 2 xzbc 
+ (a 2 — b 2 — c 2 + d 2 — 1 — a 2 — b 2 — c 2 — d 2 + 1 )y 2 + 2 xycd + 2yzbd 

+ (a 2 + b 2 - c 2 - d 2 - 1 - a 2 - b 2 - c 2 - d 2 + 1 )z 2 + 2xzbc + 2 yzbd = 0 

( b 2 + d 2 )x 2 + (b 2 + c 2 )y 2 + (c 2 + d 2 )z 2 — 2 xycd — 2 xzbc — 2yzbd = 0 
[bx — cz) 2 + (dx — cy) 2 + (by — dz) 2 = 0 
which implies that the three identities must hold: 

(6) bx = cz, dx = cy, by = dz 
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Next, identities ^ are plugged into the first equation of [ 4 ] and simplified using the 
identity a 2 + b 2 + c 2 + d 2 = 1: 

2 (a 2 + c 2 — l)x + 2 yab + 2 ycd + 2 zbc — 2 zad = 0 

2(a 2 + c 2 — l)x + 2 zad + 2 ycd + 2 zbc — 2 zad = 0 

(a 2 + c 2 — l)x + ycd + zbc = 0 

When i/O there is 

(a 2 + c 2 — l)x + ycd + zbc = 0 

( a 2 + c 2 — \)x + yc -b 2 — c = 0 

(a 2 + c 2 - l)a; 2 + y 2 c 2 + z 2 c 2 = 0 


c 2 (x 2 + y 2 + z 2 ) = x 2 (l - a 2 ) 


( 7 ) 




\/x 2 +y 2 


( 8 ) 


for a £ { — 1 , 1 }. In the other case, when x = 0 is is easy to observe from ([b]) that 

cy = 0, cz = 0 

plugging the above to Q: 

' yab — zad = 0 

(a 2 — b 2 — c 2 + d 2 — 1 )y + 2 zbd = 0 

(a 2 + b 2 — c 2 — d 2 — 1)2 + 2 ybd = 0 

a 2 +b 2 + c 2 + d 2 - 1 = 0 
replace in the second and the third equation of ([8]) with by = dz: 

( a 2 -b 2 - c 2 + d 2 - 1 )y + 2 b 2 y = 0 

(a 2 + b 2 - c 2 - d 2 - 1 )z + 2 d 2 z = 0 

(a 2 + b 2 — c 2 + d 2 — 1 )y = 0 

(a 2 + b 2 — c 2 + d 2 — 1 )z = 0 


( 9 ) 


( 10 ) 


it is impossible that both y = 0 and z = 0 since in this case x = 0 and ||u|| 7^ 0 by 
assumption. Hence, from any of the above equations there must be: 

a 2 + b 2 - c 2 + d 2 - 1 = 0 

subtracting from both sides of this equation the fourth equation of ([8]) one gets 
immediately: 

(11) (a 2 + b 2 - c 2 + d 2 - 1) - (a 2 + b 2 + c 2 + d 2 - 1) = -2c 2 = 0 

(12) =>c = 0 

which is also covered by the general solution 0 thus it can be assumed that 0 
is the only solution. Remaining spinor components are derived as follows. When 
x 7^ 0 one writes using identities (rol): 


, z 

b = —ax 


Vl^ 


VT 


sjx 1 + y 2 + z 2 sjx 1 +y 2 + z 2 


( 13 ) 
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\/l — a 2 




( 14 ) d = —ax — . = ay- , _ 

x sjx 1 + y 2 + z 2 \Jx 1 + y 2 + z 2 

in the case when x = 0 from <H3 there is c = 0 . To calculate b and d one rewrites 


( 15 ) 


' yab — zad = 0 

(a 2 — b 2 + d 2 — 1 )y + 2 zbd = 0 
(a 2 + b 2 — d 2 — 1 )z + 2 ybd = 0 
o 2 + 6 2 + d 2 - 1 = 0 


Consider the case y ^ 0 then from ([h]) there is b = Plugging this identity into 
the fourth equation of ( 15 ) one obtains: 

a 2 - 1 + (^ + l)d 2 = 0 
y z 

(a 2 - l)y 2 + (z 2 + y 2 + z 2 )d 2 = 0 


d = ay 


Vl^a? 


\/ x 2 + y 2 + z 2 


and from b = -d there is: 
v 


b = -ay- 


Vl — 


= az- 


a/1 — 1 


y \/ x 2 + y 2 + z 2 \J x 2 + y 2 + z 2 

In the case when y = 0 there must be z / 0 since ||r>|| 0 . Thus, from (Jbj) there is 

d = | b. Plugging this identity into the fourth equation of ( 151 one obtains: 

a 2 -l + (^ + l)b 2 = 0 
(a 2 - l)z 2 + {x 2 + y 2 + z 2 )b 2 = 0 


b = az 


\/T^ 


\Jx 2 + y 2 + z 2 


and from b = -d there is: 
v 


d= -az- 


V / T^2 


^~ __ = cry- _ _ 

z sj x 2 +y 2 + z 2 sjx 2 + y 2 + z 2 

In all cases a general solution to Q in coordinates is: 

( 16 ) a e [—1; 1], cr€ {-1,1} 


V / T^2 


( 17 ) b = az- 


Vl — i 


c = ax- 


Vl — i 


d = ay- 


t/1 — 


y/x 2 + y 2 + z 2 ’ y/x 2 + y 2 + z 2 ’ y \fx 2 + y 2 + z 2 

There are two parametrized solutions to the set of equations 


s =a + az 


\/l — a 2 


\J x 2 + y 2 + z 2 

vr^ 


ei2+ 


sjx 2 + y 2 + z 2 


e23 + vy 


VT^d 2 


\J x 2 + y 2 + z 2 


e3i 
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s =a + a 


s =a + a 


s =a + a 


a/ 1 — a 2 


sjx 1 + y 2 + z 2 




\/x 2 + y 2 + z 2 


{ze 12 +xe 23 + ye 3 i) 


e i2 3 {ze 3 + xei + ye 2 ) 


vT^ 


-ei 2 3 v 


s =a + aei2sV 1 — a 2 n(u) 

Finally, since many of the operations were reductions, all solutions are plugged into 
the original set of equations to check their validity. It can be seen that: 

sus -1 = (a + <rei23\/l — a 2 n(u))u(a + cre 123 \/l — a 2 n(u))^ 1 
= (a + crei 23 \/l — « 2 n(u))u(a — uei2 3 \/l - a 2 n(u)) 

= ( av + erei2 3 \/\- a 2 n(u)u)(a - crei 2 3\/l — a 2 n(u)) 

= (av + cre 123 \/l - a 2 ||u||)(a - cre 123 \/l - a 2 n(u)) 

= ava + av(—ae\ 23 \/l — a 2 n(u)) + (crei 23 \/l — a 2 ||u||)a+ 

(ae 123 \/l - a 2 ||u||)(-(jei23\/l - a 2 n(u)) 

= ava + av(—aei 23 -\/l — a 2 n(v)) + au(cre 12 3 \J\ — a 2 n(u))+ 
cr(-cr)ei23ei 23 V / l - a 2 \J 1 - a 2 ||u|| n(u) 

= ava + (1 — a 2 )v = a 2 v + v — a 2 v = u 

which conhrms that all the solutions are valid. Now we obtain the following formula: 

s = a + ae 123 y/l — a 2 n(v) 

for [—1; 1]. Now it is further simplified by using trigonometric and exponential 
series converging for all arguments. We substitute a = cos(a) to obtain: 

a + ae 123 \/l — a 2 n(u) = cos(a) + sin(a)ei 23 n(u) 

where a £ [ 0 ; 27 r). Then, we express the trigonometry using exponential function. 
We use series expansions which are valid for all a. Note that (—l) fc = (e 12 3e 12 3) fe = 
(e 123 ) 2k and (n(u)) 2fe = 1. 


k , „ ( -1 ) fc „.2fc + l 


cos(a) + sin(a)e 12 3 n(u) = ^ + ei 23 n(«) ^ (2fc + 1) , 

(ei23) 2fc 


a 


k =o 




\2k 


k =0 
oo 


(2fe)! 




= £ 

k—0 

oo 

= £ 


\2 k 


oo / \2k-\-l 

V e 123J 2fc+l / / \2fc+l 

—„ („M) +y. (2t+1)! a <"(>-» 


( e 123) 2 k 


(e 123 an(v)) 2k + ^ (ei 23 an(u )) 2fc+1 


fc=o 


(2fc)! 


fc=0 


(2/c + l)! 



SWING-TWIST DECOMPOSITION IN CLIFFORD ALGEBRA 


9 


( e i23an(i;)) fc 

= 2^- J 7 \ - = exp(ei 23 a n(i/)) 

k =0 

which is precisely the stated general formula. This completes the proof. □ 

5. A SPINOR WHICH ROTATES A GIVEN VECTOR 

The following lemma states what is a spinor which rotates a given vector to be 
equal to a target one. In literature it is sometimes called direct rotation , as in 
[ 3 ]. Schematic view is presented in figure [3j The lemma is required to prove a 



Figure 3 . A direct rotation of AB to AB' 

general theorem about the set of spinors which rotate a given initial vector to a 
given target vector. We start with the following initial lemma which applies to the 
Clifford algebra of R 2 . At first the following basic property of Cl 2 is proved: 

Lemma 1 . The external product of two given vectors v and w in Ci 2 can be written 
as: 

w Av = - sgn(w A wei2)||r> — u>|| ||u + u>||ei2 

Proof. Assume that v = a;ei +ye2 and w = pe i + ge2. In special case when w = Av 
both sides of the equation are equal to zero therefore the equation holds. In general 
case w ^ Av and one can rewrite both sides in coordinates. Left-hand side is equal 
to: 

w A v = -(wv - vw) = - ((pei + qe 2 )(xe 1 + ye 2 ) - (xe ± A ye 2 )(pe i A qe 2 )) 

= ^ ( 2 pye 12 - 2 qxe 12 ) = (py - qx)e 12 
Right-hand side is equal to: 

- sgn(i; A wei 2 )— HHI' 1 ’ + R’|| e i2 
= ^ sgn(^(i;u) - wv)e 12 )y/(v - w) 2 ^/(v A w) 2 e 12 

= ^ sgn(((a;ei + ye 2 )(pe 1 A qe 2 ) 

- (pe i +ge 2 )(a;ei + pe 2 ))ei 2 ) yj(v - w) 2 (v A w) 2 e 12 
= ^ sgn((2pa;ei2 - 2 pye 12 )e 12 ) 

V(( x - P ) 2 + (y ~ q) 2 )((x A p ) 2 A(yA q) 2 )e 12 
= ^ sgn (py - qx)yjA(py - qx) 2 e 12 = sgn (py - qx)jpy - qx\e X2 
= (py ~ <ix)ei 2 

Both sides are equal so it completes the proof. □ 
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Lemma 2. Let v and w be a pair of non-zero vectors of equal length in C ^2 such 
that v ^ —w. Then there exist a pair of spinors ±s £ Spin( 2 ) which rotate v so 
that the rotated vector is equal to w: 

svs ^ 1 = w 

then the pair of spinors is equal to: 

s = ± n(u + w) n('u) 

Proof. Let s = a + be 12 £ Spin( 2 ) be a spinor and v = xei + ye 2, w = pei + qe 2 
be a pair of given vectors. The equation can be expressed in coordinates with the 
following set of equations: 

f ( a 2 — b 2 )x + 2 aby = p 
1 (a 2 — b 2 )y — 2 abx = q 


( 18 ) 


From the assumption it is impossible that both x and y are simultaneously zero. 
When i/Owe solve the set of equations: 


2 2 P~ 2ab V 

a — b = - 

x 

p — 2 aby 

- y — Zabx = q 

x 

The second equation is then simplified: 

(p — 2 aby)y — 2 abx 2 — qx = 0 
py - qx 


2 ab = 


2T 


Which is then plugged into the first equation of ( 18 ): 


(a 2 - b 2 )x + 


py - qx 


(2a 2 — l)x = p — 


2 V=P 

-y 2 

py 2 - qxy 


x* 


r 


2 _ px 2 + py ' 2 - py 2 + qxy 


2a 2 = 


a = ±1 


x(x 2 + y 2 ) 


+ 1 = p_E±M + 1 


I x 2 + y 2 + px + qy 


V x*+y* 

The same result can be obtained when we consider the other case y 7^ 0 : 

q + 2 abx 


2 l2 

a — b = 
q + 2 abx 


y 


y 


x + 2 aby = p 


The second equation can be simplified: 

(q + 2 abx)x + 2 aby 2 — py = 0 
py - qx 


2 ab = 


r 
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Which is then plugged into the second equation of ( 181 : 


(a 2 -b 2 )y- 


py-qx 
x 2 + y- 


; x = q 


pxy — qx 
x 2 + y 2 

2 qx 2 + qy 2 + pxy — qx 2 


( 2 a 2 - 1 )y = q + 
2 a 2 = 


y(x 2 + y 2 


+ l = E|±w + l 

ar + y z 


a = ±4 


I x 2 + y 2 + px + qy 
2 (x 2 + y 2 ) 


We obtained one general solution which is correct for all x, y. Now, for each 
a € {—1,1} of a = ~2(x^+y^) qV we h &ve exactly one corresponding solution 

b. It can be calculated by plugging it into 2 ab = 2 (x^+y 2 ) v& lid f° r a ll x - V■ When 
b/0 then: 


b = 

b = 

b = a 


b = a 

b = a 

b = a 

b = a 


py - qx 

2 a(x 2 + y 2 ) 

(py - qx)\j 2 .(x 2 + y 2 ) 

2 a\J x 2 + y 2 + px + qy(x 2 + y 2 ) 
_ PV ~ qx _ 

yjx 2 + y 2 +px + qy^j 2 (x 2 + y 2 ) 
(py - qx)sjx 2 + y 2 - px - qy 


sj(x 2 + y 2 ) 2 - (px + qy) 2 ^j 2 (x 2 + y 2 ) 

(PV ~ qx)\/x 2 +y 2 -px-qy 
\J(x 2 + y 2 ) 2 - (px + qy) 2 \j 2 (x 2 + y 2 ) 

_ (py - qx) yjx 2 + y 2 - px - qy _ 

\]x^ + y A + 2 x 2 y 2 — p 2 x 2 — q 2 y 2 — 2 pqxy^j 2 (x 2 + y 2 ) 


(py ~ qx)\/x 2 + y 2 — px — 


qy 


\Jx i + y 4 4- 2 x 2 y 2 — (x 2 + y 2 — q 2 )x 2 — (x 2 + y 2 — p 2 )y 2 — 2 pqxy 

1 


\J 2 (x 2 + y 2 ) 

py - qx 


b = a 


I x 2 + y 2 — px — qy 


A /(py - qx ) 2 V 2 (x 2 + y 2 ) 


note that (py — qx)e \ 2 = wAv but it is assumed that v is not parallel to w therefore 
py — qx 0 . Since e 2 2 = —1 the following applies: 


py ~ qx 
a/ (py - qx¥ 


= sgn(n A w;ei2) 


and b is equal to: 


b = a sgn(i; A uie 12 ) 4 


I x 2 + y 2 — px — qy 
2 (x 2 + y 2 ) 
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When a = 0 then from b 2 = 1 — a 2 we know that b = ±1 which is also a special 
case of the general solution. Finally, we simplify the formula to a coordinate-free 
solution. For each a £ { — 1 , 1 } there is a solution: 

x 2 + y 2 + px + qy . . x 2 + y 2 — px — qy 

8 = "V 2(i 2 + j/ 2 ) +aSgn{V A " ,eia) V 2(a: 2 + y 2 ) ^ 

8= y2 ( J + , 2 ) (v,l2+1,2+pl+ro+ 

sgn(u A uiei 2 ) \/ x 2 + y 2 — px ~ qye i2 ) 


V / 2( J+ y2) ( V 5 (x2 + y 2 + P 2 + g 2 + + gy)) + 


sgn(u A we 12 )\l - (x 2 + y 2 + p 2 + q 2 - 2(px + qy))e 12 ) 


s = 


2yjx 2 + y 2 


{y/x 2 + y 2 +p 2 + q 2 + 2(px + qy)+ 


sgn(v A wei 2 )yjx 2 + y 2 + p 2 + q 2 — 2(px + qy)e i2 ) 


s = 


2 v 


-{y/w + WW + VW + WV+ 


sgn(u A wei 2 )y/vv + ww — vw — wve 12 ) 
s = 


2||u 

(J 


2\\v 


{yj(v + w) 2 + sgn(u A we r2 )y/(v - w) 2 e 12 ) 
(\\v + w II + sgn(u A we 12 )||u - w||e 12 ) 


From the assumptions, there is ||u + w || ^ 0 so: 


s = 


\v + HI 2 + sgn(u A ruei 2 )||u — HI ||u + ic||ei 2 ) 


2||u|| ||u + HI 

Next, the basic identity from lemma[T]is used so the formula can be rewritten with 
only external product: 

w A v = — sgn(u A wei 2 )||u — Hl||v + Hl e i 2 

It reads that the external product is equal to a bivector whose area is equal to the 
half of the area of the parallelogram spanned by the both vectors and with the sign 
adjusted to the sign of the external product. Using the above identity, the following 


(\\v + ru|| 2 + 2v A w) 

(vv + ww + 2v ■ w + 2v A w) 


s = 


s = 


are 

possible: 

a 

2\\ V 

||v + it; || 
a 

2\\v 

II ll' y + HI 


(■ 

INI 

||u + w || 




||u|| \\v + it;11 
s = ± n(u + w) n(v) 
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Which is precisely the stated formula. □ 

We observe that the same vector formula applies to any Clifford algebra of n- 
dimensional real space. Thus, we propose the following: 

Proposition 2. Let v and w be a pair of non-zero vectors of equal length in C t n 
such that v —w. Then there exist a pair of spinors each of which rotate v so that 
it is equal to w: 

svs~ l = w 


the pair of spinors is equal to: 

s = ± n(i> + w) n(i>) 

Proof. Instead of the original equation svs^ 1 = w an equivalent equation will be 
proved: 

sv = ws 

Let v and w are the given vectors. Start with: 

v — v = v — v 

since vv = ww > 0 each element can be multiplied by the scalar vv or ww: 


vvv — wwv = ww — ww 
vvv + ww = ww + wwv 
(v + w)vv = w(v + w)v 


it is assumed that v y^ — w so ||n + w | > 0 and both sides can be divided by 
||n + to|| \\v\\: 

n(u + w) n(v)v = w n(v + w) n(u) 

To ensure that s is a spinor it is sufficient to check whether it is a direct product 
of a scalar and a bivector and whether its norm is equal to 1 . Indeed: 


s = n(w + w) n(u) 


vv + wv 
||i) + w|| |M| 


||t>|| 2 + wv 
||i> + u>|| ||u|| 


G /\ M 3 (g) /\ M 3 
0 2 


and since 


s 1 = n(||u|| 2 + w ■ v + w A v) 1 

= n(||u|| 2 + v ■ w + v Aw) = n(vv + vw) = n(t>) n(i> + w) 

the norm is: 

ss 1 = n(v + w) n(u) n(u) n(u + w) = || n(i; + it?)|| 2 11 n(f)|| 2 = 1 
Which completes the proof. □ 


Since the theorem is valid for Cl n it is also valid for C £3. In more general 
conclusion, for any dimension there exists a simple formula giving a pair of spinors 
rotating a given vector so that it covers another given vector of the same length. 
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6 . The set of spinors which rotate a given initial vector to a given 

TARGET VECTOR 

There are two ways of defining a swing-twist representation (composition or 
decomposition): 

• twist q is done before swing p; for a given spinor r we have: s = pq. 

This is the swing-after-twist representation. 

• twist q is done after swing p; for a given spinor r we have: s = qp. 

This is the twist-after-swing representation. 

The swing-after-twist representation is used in [T] while the twist-after-swing rep¬ 
resentation is usually used in the direct method. Since in practice both conventions 
are used (the first one is slightly less common) in this paper both decompositions 
will be presented and proved. The following theorem is a general solution to Clifford 
product equation: 

si>s -1 = w 

We prove the following 

Proposition 3. Let v and w be a pair of non-zero vectors of equal length in CC 3 
such that v 7^ —w. The set of spinor solutions s £ Spin( 3 ) to 

sus -1 = w 

in the case of swing-after-twist representation is equal to: 

s = ± n(u + w) n(u) exp(ei23a n(u)) 
and in the case of twist-after-swing representation is equal to: 

s = ± exp(e 12 3a n(u>)) n(u + w) n(u) 

Proof. The solution is the set of all possible rotations which swing axis from initial 
v to terminal w with any possible twist during the movement. Using propositions [l] 
and [2] one composes swing and twist according to the order used in a given represen¬ 
tation. In the case of swing-after-twist representation twist q = exp(e 123 o: n(w)) fac¬ 
tor (around v axis) is applied at first and only after it, swing factor p = n^c+w;) n(u) 
is applied (swinging the axis from v to w). In the other case of twist-after-swing 
representation swing p = n(v + w) n(u) is applied first (swinging the axis from v to 
w) and after that twist factor q = exp(ei 2 3a n(ic)) is applied with respect to the 
terminal axis w. □ 

7. Swing-twist decomposition of a spinor 

This is the main result of this paper. In this section the inverse of formulas given 
in proposition [ 3 ] is derived. Given a spinor it is possible to calculate its decomposi¬ 
tion into twist and swing factors in respect to a given non-zero vector. Since there 
are two different swing-twist representations, there are also two swing-twist decom¬ 
positions for swing-after-twist and twist-after-swing representation respectively. In 
this section, the initial vector will be called a base vector. 

Theorem 1 (Swing-twist decomposition of a spinor in swing-after-twist represen¬ 
tation). Assume that s £ Spin( 3 ) is a spinor. For any non-zero base vector v £ C^3 
such that svs^ 1 ^ —v there exist a unique up to the sign swing-twist decompo¬ 
sition in swing-after-twist representation 

s = ±pq 
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where swing spinor p and twist spinor q are equal to: 

p = ±sa v (s) 
q = ±<r v (s) 

where a v (s) : Spin( 3 ) —> Spin( 3 ) is a function of spinor s: 

<Jb{s) = n(v(v • s)) 

the reversion a v (s ) is equal to: 

cr v (s) = n(v(v • s)) 

The function a v (s) will be called a twist projection function. 

Proof. The proof is divided into several steps. At first the problem is reformulated 
in coordinates. Then, twist angle is calculated and twist spinor and finally swing 
spinor. Denote the following: 

s = a + be 12 + ce 2 3 + cie 3 i 
v = Ua-ei + v y e 2 + v z e 3 
w = u^ei + w y e 2 + w z e 3 
The angle a will be described implicitly by: 

cos(a) = k 
sin(a) = l 
k 2 + l 2 = 1 

Vector w exists on the assumption that sus -1 ^ —v. What does this assumption 
require about the base vector and the spinor is explained in the appendix of this 
paper. At first, the equation is rewritten in coordinates. The swing factor is equal 
to: 


p = n (w + v) n(u) 

_ ((w x + Ut)ei + ( w y + v y )e 2 + (w z + n z )e 3 )(n a ,ei + v y e 2 + v z e 3 ) 
\/(w x +v x ) 2 + (Wy +Vy) 2 + (w z +v z ) 2 ^vl + U.2 + v\ 

= i(w x + V X ) 2 + (Wy + v y )' 2 + ( w z +v z ) 2 )~i{vl + Vy+vl)~i- 
[K + v x) v x + ( w v + v y) v v + i w z + V z )v z + 

((w x + V X )Vy - (Wy + Vy)v x )e 12 + ((Wy + Vy)V Z ~ (W Z + V Z )Vy)e 23 + 

((w z + v z )v x - ( w x + v x )v z )e 3 i] 


Twist factor is equal to: 

( 19 ) q = cos(a) + e 123 n(u) sin(a) = k + 



(b z e 12 + b x e 23 + 6 y e 31 ) 


Combining the above formulas there is: 

pq = (K +v x ) 2 + (w y + v y ) 2 + {w z + v z ) 2 )~^ {v 2 x +v 2 y +v z )~5- 

[( W x + V X )V X + (Wy + Vy)Vy + (W Z + V Z )V,+ 

((w x + V X )Vy - (Wy + Vy)V X )e 12 + ((Wy + Vy ) V Z - (W Z + V Z )Vy) e 23 + 

((w z + v z )v x - (w x + v x )v z )e 3 i}- 
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[fc- 


y 


V 2 + V 2 + V 2 


~( b z e 12 + b x e 23 + &. y e3i)] 


= (( w x + v x ) 2 + (Wy + Vy ) 2 + ( w z + v z ) 2 ) 2 ( vl + Vy + Vz ) 2. 
+ V x ) v x + ( Wy + Vy)Vy + (W Z + V Z ) V Z )+ 


( k ( 


W X Vy ~ W„V : 


y 


) + l\Jvl + o 2 + v 2 ,{w z + v z ))e 12 + 

( k(w y v z - w z v y ) + lyjv% + v 2 + v 2 (w x + v x ))e 23 + 

(k(w z v x - w x v z ) + l\Jvl + v 2 + v 2 (w y + v y ))e 31 ] 

Two spinors are equal if and only if corresponding coefficients are equal. There¬ 
fore the following set of equations determines the solution: 

'a = [ k((w x + v x )v x + (w y + Vy)v v + (w z + v z )v z )]A 
b = [k(w x Vy - w y v x ) + lyjv* + v 2 +v^(w z + o z )]A 
c = [k(w y v z - w z Vy) + l^jv 2 + v 2 + v 2 (w x + o x )]A 
d = [ k(w z v x - w x v z ) + lyjv% + o 2 + v^(w y + v y )\ A 
fc 2 + l 2 = l 


( 20 ) 


w 2 x + w 2 y + to 2 = vl + v 2 + v 2 


[ a 2 + b 2 + c 2 + d 2 = 1 


Where A = ((w x +v x ) 2 + (w y +v y ) 2 + (w z +v z ) 2 )~2 (f x +u 2 +u 2 ) _ 5 . By substituting 
variables in the set of equations (201 it is easy to rise unfavourably the degree of 

ving careful 
k from the i 

k(w + v ) ■ v 


involved polynomials. Thus the following careful operations are performed. First, 
note that it is easy to obtain variable k from the first equation: 


( 21 ) 

Observe that: 


\w + u v 


1 + u|| 2 = (w + v)(w + v ) = ww + vv + wv + vw = 2 (v • v + v ■ w ) = 2 v ■ (w + v ) 


Using this, one simplifies ( 21 ) to: 

k(w + v ) ■ v 
||to + r'|||H| 

From that one obtains the value of k : 


( 22 ) 


k = 


k\\w + v \\ 2 _ \\w + u| 

|to + v|||M ' 211-^11 


2|h 


to + v 


The value of k is now put into the set of equations ( 20 ) which then can be simplified 
to: 


( 23 ) 


' b\\w + a|| = 2 a ||t|) ^n ( w x v y - w v v x ) + l ( w z + v z ) 
c||to + o|| = 2a ( w y v z - w z v y ) + l { w x + v x ) 
^ d\\w + 'o|| = 2a ,, * ( w z v x - w x v z ) + l ( w y + v y ) 
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It is a set of three linear equations with three unknowns M[w x , w y , w z ) T = N. The 
characteristic matrix M is: 


M = 


2a i 


— 2a 


'IU+'U|| 

l 

l 


—2a | 
2a-rr 


Ul+'l'll 


H-ull 1 

l 


—2a j 
2a i 


u’+i'll 




and the vector N is equal to: 


N = [b\\w + v|| - lv ~, c||ui + u|| - lv x , c?||w + u|| - lv y ] T 


Any known method can be used to compute the determinant of M which is equal 
to: 

det(M) = l 


The case l = 0 needs special care and will be now discussed separately. In this case, 
from ( 19 ) twist factor is equal to: 

l 


q = k + 


-~{b z e 12 + b x e 23 + &j,e3i) — k 


Since q is a spinor, its component k must be equal to: 


k = ± 1 


which then implies that twist and swing factors are equal to: 

p = ±s 

q = ± 1 


It is a special case of the general formula. It holds for all vector v and spinors s 
such that a = zir, z £ Z. In the case of l 7^ 0 there exist exactly one solution to 
( 23 ). After having the matrix inverted and solution calculated, one gets: 


Wx = [ 4a 2 (vt+vl+vt)l + \\w + vfl 3 }- 1 

[2a(bv y - dv z )\\w + u|| 2 Z+ 

||w + u|| 2 Z 2 (c||'u> + v\\ - v x l) + 4 a 2 v x ((bv z + v x c + v y d)\\w + u|| - (v 2 x + v 2 + v 2 z )l)\ 
w y = [ 4 a 2 (u 2 +v y + v 2 z )l + \\w + u|| 2 Z 3 ] -1 
[2a(— bv x + cv z )\\w + v\\ 2 l+ 

||w; + u|| 2 Z 2 (d||ui + u|| - v y l) + 4 a 2 v y ((bv z + v x c + v y d)\\w + u|| - (u 2 + v 2 + v 2 z )l)\ 
Wz = [ 4 a 2 (v 2 x + v 2 y + v 2 )l+\\w + v\\ 2 l 3 ]- 1 
[2a(cv y + dv x )\\w + u|| 2 Z+ 

\\w + v\\ 2 l 2 {b\\w + u|| - v z l) +4 a 2 v z (( bv z + v x c + v y d)\\w + v\\ - (v 2 + v 2 + v 2 z )l)] 


All three denominators are equal and can be simplified: 


4 a 2 (u 2 + v 2 + v 2 )l + ||w + u|| 2 Z 3 

= Z( 4 a 2 |M| 2 + ||w + u|| 2 (l - a 2 /y )) = l\\w + u|| 2 

||ui + v\\ 2 
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A compact form is achieved with the following simplifications: 


w. 


1 


I w + v\\ 2 


(2 a 


bv y — dv z 
cv z - bv x 

dv x — CVy 


c||w + i'll - v x l 
|u' + v\\ 2 l + ||«) + i'll 2 / 2 (i||w + w|| - v y l 
b\\w + i>|| - v z l 


4 a 2 


(( bv z + cv x + dv y )\\w + v\\ - ||p|| 2 0 ) 



V x 


c 


c 

2 a 

Vy 

X 

d 

+ l\\w + V I 

d 


v z 


b 


b 


- I 


4 a 2 


|w + d | 2 


V x 


Vx 


c 

Vy 

(lk + ^ll 

Vy 


d 

Vz 


V z 


b 


- V 


Vector [c, d , b] T can be formulated with coefficients of spinor s: 

*[s]2 = — ei 23 [s ]2 = —6123(^12 + ce 23 + de3i) = cei + de 2 + be 3 
With this substitution further simplification is possible: 

4 a 2 


w = 2av x *[s] 2 + l\\w + v || * [s]2 — l 2 v + 
but form previous equations there is: 

( 24 ) ;2 = 1_a2 I 

so after substitution one obtains: 


Z||w +1;|| 2 
4 |H| 2 

w + v || 2 


v(\\w + V\\V ■ * 8 2 - u 


= 2 av x *[s] 2 + l\\w + u|| * [a]2 — v + 


4 a 2 v 
|w + v\ 


;v-*\s 2 


The above formula is useful to compute ||ie +i>||. Note that one can move — v to 
the left hand side of the equation and compute lengths of both sides. 


\w + u|| 2 = (2 av x *[s] 2 + l\\w + v || * [s] 2 + 


4 a 2 v 


\w + i'll 2 


v ■ *[s] 2 ) 2 


4a 2 v 


= (2 av x *[s] 2 ) 2 + (l\\w + V || * | S ] 2 ) 2 + ( ... ^ N2 i> • *[s] 2 ) 2 

l\\w + or 


+ 2[(2 av x *[s] 2 ) • {l\\w + u|| * [s] 2 ) + (2 av x *[s] 2 ) • 


4 a 2 v 


\w + u || 2 


V ■ *[s] 2 ) 


+ (!||io + v|| * [«]j) • (- *[*W1 


= 4 «> X *[»] 2 ) 2 + ! 2 »«> + » f (*[» l 2) 2 + ■ *[»] 2 ) 2 + 8 a > • *[»] 2 ) 2 

l Z \\W + v\\ z 

From the property that (&x*[s]2) 2 + (6-*[s]2) 2 = IH| 2 (*[s]2) 2 the following formula 
is obtained: 

|| W +u|| 2 =4a 2 ||b|| 2 H s ]2) 2 +/ 2 || W +u|| 2 (x[ s ] 2 ) 2 +4a 2 (u-x[ s ] 2 ) 2 + 

l z \\w + v\\ z 
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There is: 

( 25 ) l' 2 \\w + v\\ 2 = fl - * a f 1 ' ) |k + -f = |k + v\\ 2 — 4 a 2 ||t>|| 2 

V \\w + v\\ 2 J 

So one can write the following equivalent formula: 

Ik + = 4 a 2 |H| 2 (*[ S ] 2 ) 2 + (||u, + nil 2 - 4 a 2 |M| 2 )(*[ S ] 2 ) 2 + 

. 9 . r i \2 16 a 4 || vll 2 , r i \2 

4 a (v*[fl] 2 ) + Z2 || w + ,'||2 (--*M2) 

= \\w + u || 2 (*[ S ] 2 ) 2 + 4 a 2 (i; • *[s ] 2 ) 2 + • */] 2 ) 2 

r\\w + v\\ z 

Substituting (*[s ] 2 ) 2 = 1 — a 2 and once again | 25 | ) one obtains: 

||ti) + v\\ 2 = ||ui + n|| 2 (l - a 2 ) + (v ■ w[s ] 2 ) 2 
Which can be simplified to: 


4 a 2 (||u; + u || 2 — 4 a 2 ||n|| 2 ) + 16 a 4 ||i ;|| 2 
l 2 \\w + i >|| 2 


a 2 ||w + u || 2 = {v -*[s] 2 ) 2 ^- 
a 2 (||w + v\\ 2 - 4a 2 ||u|| 2 ) = 4 a 2 {v ■ *[s ] 2 ) 2 
There are two cases. When a / 0 the following equation applies: 

Ik + -l| 2 = 4( , c • *[s] 2 ) 2 + 4a. 2 |k|| 2 
Using the above identity, there is: 

l 2 \\w + u|| 2 = ||tu + z;|| 2 - 4a. 2 11^11 2 = 4(v • w[s] 2 ) 2 + 4a 2 ||u|| 2 - 4a 2 ||/| 2 
Z 2 |k + u || 2 = A(y ■ *[s ] 2 ) 2 
Modifying the identity p4| there is: 


Z 2 = 1 — a 2 - 


4|h 


Ik + -II 2 

|tu + v|| 2 — 4a 2 ||u|| 2 = Z 2 ||w + t>|| 2 
k + -l| 2 — 4a 2 ||u|| 2 = 4(t> • *[s] 2 ) 2 

Ik + -II 2 = 4(v • a[s] 2 ) 2 + 4a 2 ||i>|| 2 
Ik + -II = 2a /{v ■ *[s] 2 ) 2 + a 2 ||v|| 2 


Using the formula ( 22 ) one obtains the value of k: 

2II-II „ 


( 26 ) 


k = i 


2|h 


Ik + -II 2 v /(io*[ S ] 2 ) 2 + a 2 |M| 2 v 2 (- - *[ s ] 2 ) 2 + a 2 IM| 2 


Again from ( 24 ) it is possible to calculate l: 

4 a 2 ||i >|| 2 


P = 1 - 0 » 4 <HI 2 =1 _ 

Ik + -II 2 4 [(u • *[s] 2 ) 2 + a 2 ||i>|| 2 ] 
(v • *[s] 2 ) 2 + a 2 


(- • *W 2 ) 5 


(- • *[ s ] 2 ) 2 + a 2 |/|| 2 (v ■ *[s ] 2 ) 2 + a 2 ||u || 2 (v ■ *[s ] 2 ) 2 + a 2 ||i >|| 2 

It is possible to take the square root and abandon sign. This comes from the fact 
that when s is a given spinor then the solution represented by variable Z must 
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include the solution for — s as well. In that case the numerator is negated, so also 
—I is a solution automatically. 


( 27 ) 


l = 


v • *[s ] 2 


vV*H 2 ) 2 + a 2 HI 


In the case when a = 0 from ( 24 ) there is: 

/ 2 = 1 

k 2 = 1 - l 2 = 0 


The same argumentation as above can be used to select any solution of l 2 = 1 . 
Both cases are special cases of ( 26 ) and ( 27 ) so it is sufficient to consider only those 
general solutions. The twist factor equal to: 


q = k + ei23 n(v)l = 

= _aJHj_ e 123 n(t>)(p • *[s] 2 ) 

\J( v ' *[ s ] 2 ) 2 + 4a 2 ||w|| 2 \J(v ■ *[s] 2 ) 2 + a 2 |H| 2 
a||u|| + n(ri)(ri • *[s]) _ a||v|| + n(v)(v ■ s ) — n(n)(n ■ [s] 0 ) 

“ V( v ■ *W 2 ) 2 + _ a 2 lR 2 ~~ V(v • *[s] 2 ) 2 + a 2 lRF 

a||u|| + n(u)(v • s) — a|MI n(p)(p • s) 

= y/(v ■ *[a] 2 ) 2 +^RF = TCT^F+^RF 
Define a norm of a pinor p = p 0 + pi 2 ei2 + P23&23 + ^31631 G Pin( 3 ): 


n(p) := 


P 

Vpp 


Po + Pl 2 ei2 + P2 3 e 23 + P31 e 31 
VPo +Pl2 +P23 +P31 


€ Spin( 3 ) 


where p = p 0 - pi 2 ei 2 - p 23 e 23 - p 31 e 31 . 

Using the fact that s = a+ei 2 3*[s ] 2 the numerator of twist factor can be written 
as: 


n(v)(w • s) = n(w)('y • (a + e 123 * [s] 2 ) = n(w)(z;a + e 123 v • *[s] 2 ) 

( 28 ) = a|H| +ei 23 n(w)t; • *[s ] 2 

the norm of the numerator of twist factor is equal to: 

V HMI + e i 23 n(n)p • *[s] 2 )(a||t;|| - e 123 n(u)p • *[s] 2 ) 

= V « 2 ||^|| 2 + (ei 23 n(v)v ■ *[s ] 2 ) 2 = Va 2 \\v\\ 2 + (v • *[s ] 2 ) 2 

one notices that it is equal to the denominator. In result, it is possible to write 
simply: 

q = n(v(v ■ s)) 

The swing factor is calculated by modifying the original equation: 


s=pq 
p = sq~ 


-1 


by using (281 twist inverse can be easily calculated: 

q ^ 1 =q = a v (s) = a|H| + ei 23 n(v)v ■ *[s ] 2 = a||w|| - ei 23 n(v)v ■ *[s ] 2 
= a||f || + e i23 n(ri)n • *[-s ] 2 = n(v)(va + ei 23 u • *[-s] 2 ) 

= n(v)(v ■ (a + e 123 * [-s] 2 )) = n(v)(v ■ s) 
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where a v (s) = n(w)(t; • s). Both swing and twist factors can be negated as a given 
spinor and its negation define the same rotation: 

p = ±stT„(s) 
q = ±cr v (s) 

which completes the proof. □ 

The following dual theorem can be easily proved with the previous theorem 0 : 

Theorem 2 (Swing-twist decomposition of a spinor in twist-after-swing represen¬ 
tation). Assume that s £ Spin( 3 ) is a spinor. For any non-zero base vector v £ C^3 
such that svs 7^ — v there exist a unique up to the sign swing-twist decompo¬ 
sition in twist-after-swing representation 

s = ±qp 

where swing spinor p and twist spinor q are equal to: 

p = ±o v (s)s 
q = ±cr v {s) 

where cr v (s ) is a twist projection function. 

Proof. Assume that u = s -1 . According to theorem [2] there exists a swing-twist 
decomposition of spinor u in swing-after-twist representation: 

u = ±pq 

where p is a swing factor and q is a twist factor of u in respect to the base vector 
v. Taking inverse of both sides one obtains: 

ir 1 = s = ±q~ l p~ l 

which is a twist-after-swing decomposition of spinor s. Swing and twist factors can 
be rewritten as: 

p -1 = = (S - 1 CT„(S -1 )) -1 = <T„(S -1 ) -1 S = CT„(s)s 

q — cr v (vj — (Tyi^s ) o’y(s') 

which completes the proof. □ 

8. Twist projection function 

In this section basic properties of the twist projection function are discussed. 
Intuitively, twist projection function takes a spinor and returns its twist factor in 
respect to a given vector. In this paper the following definition was assumed 

a v (s) := n(v(v ■ s)) 

This is algebraically equivalent to projecting it onto a vector in the sense of Clif¬ 
ford algebra. This an improvement over the work of Huyghe [Tj where the author 
introduces similar projection operator artificially in quaternion algebra. Twist pro¬ 
jection function is a projection due to the following 

Proposition 4. A twist projection function is a projection. 
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Proof. It is enough to prove that <j v (a v (s )) = a v (s) for any vector v and any spinor 
w: 


o- v (cr v (s )) = n(u(t> • n(v(v • s)))) = n(u(u 


a|MI + e 123 n(ll)(l> ■ *[s] 2 ) 
\J( v ' *[s ] 2 ) 2 + a 2 ||u || 2 


)) 


_ , V- a||ti|| + ei 2 3 (v • n(v))(v • *[s] 2 ). _ ,v- a||i>|| + e i2 3 ||'^||v ■ *[s] 2 . 

' n r t \o i on 11 o ' ' n r i \o i oTi 11 o ' 


= mv- 


V>-*[ S ] 2 ) 2 +a 2 |M | 2 

|w||v)-(a + *[a] 2 ) =)=n( || t) ||a 


\/(v ■ *[s] 2 ) 2 + a 2 ||v|| 2 
n(u)(v • s) 


VV*[ S ] 2 ) 2 + a 2 |M| 

= n(u(v • s)) = a v (s) 


V( v ' *i s h) 2 + a2 \\v\\ 


f) 


which completes the proof. 


□ 


9. Applications 

In this section an exemplary application of the proposed theoretical results is 
presented. It is a fast, concise and numerically stable algorithm for calculating 
swing-twist decomposition of a spinor. 

9 . 1 . An efficient algorithm for swing-twist decomposition of a spinor. To 

propose an efficient method of calculating swing-twist decomposition the following 
formula is used for twist projection function: 

a v (s) = n(v(v ■ s)) = n(v(v ■ (a + e 123 * [s] 2 )) = 

= n (v(va + e 123 v • *[s] 2 )) = n(a||v || 2 + e 123 v(v • *[s] 2 )) 

If a rotation is represented by a pinor normalization can be omitted. In case of 
spinor representation normalization is needed and requires square root computa¬ 
tion or equivalently requires arithmetic with square root extension. The proposed 
method uses swing-after-twist representation and is presented in algorithm |Tj 


Algorithm 1 Compute swing-twist decomposition of a spinor s = 
Require: v = xe 3 + ye 2 + ze 3 , s = a + be 12 + ce 23 + de, 

Ld -vi mi v* ^~\ • I I n i I I —/— Q 


pq 


Ensure: ||i> 

u 4— xc + yd + zb 
n t— x 2 + y 2 + z" 
mean 
l <— \frn- 

q 


•31 


m 

T _r " 
p «— sq 

return p 


+ u 2 n 

zu 0 

X e 12 - 


xu _ I yu 

T e 23 + X e 31 


Because of its simplicity the proposed method favourably compares to existing 
methods for computing swing-twist decomposition (as presented in section [2j 

Appendix A. Which decompositions are impossible 

When a combination of a spinor and a base vector cannot be decomposed into 
swing and twist. From theorem [I] or [ 2 ] a decomposition is impossible when: 

svs~ l = —v 
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the above condition can be rewritten in terms of spinor coordinates. Using © one 
writes 

' (a 2 — b 2 + c 2 — d 2 + l)x + 2 y(ab + cd) + 2 z(bc — ad) = 0 

(a 2 — b 2 — c 2 + d 2 + 1 )y + 2 x{cd — ab) + 2 z(bd + ac) = 0 

(a 2 + b 2 — c 2 — d 2 + 1 )z + 2 x(bc + ad) + 2 y{bd — ac) = 0 

a 2 + 6 2 + c 2 + d 2 - 1 = 0 


(29) 


Denote as TU, the left side of All equation of 29 A valid solution (a, b, c, d) must 
satisfy the following equation: 

xW\ + yW 2 + zWj, + ( x 2 + y 2 + 2 2 )Vb 4 = 0 
the above formula can be expanded: 

(a 2 — b 2 + c 2 — d 2 + l)a; 2 + 2 xy(ab + cd) + 2 xz{bc — ad) 

+ (a 2 — b 2 — c 2 + d 2 + l)?/ 2 + 2 xy{cd — ab) + 2 yz(bd + ac) 

+ (a 2 + b 2 — c 2 — d 2 + l)z 2 + 2 xz(bc + aci) + 2 yz(bd — ac) 

+ ( x 2 + y 2 + z 2 )(a 2 +b 2 +c 2 + d 2 - 1) = 0 

(u — 6 c — d~ U 1 A u y -f c -f- d“ — l)x~ -f- 2 xycd -t- 2 xzbc 
+ (a 2 — b 2 — c 2 + d 2 + 1 + a 2 + b 2 + c 2 + d 2 — 1 )y 2 + 2 xycd + 2 yzbd 

+ (a 2 + b 2 - c 2 - d 2 + 1 + a 2 + b 2 + c 2 + d 2 - l)^ 2 + 2 xzbc + 2 yzbd = 0 

(cx + dy + bz ) 2 + a 2 (x 2 + y 2 + z 2 ) = 0 
which implies that the two identities must hold: 

a 2 (x 2 + y 2 + z 2 ) A cx + dy + bz = 0 
but since ||u|| ^ 0 it must be 

a = 0 A v ■ w[s] 2 = 0 

These conditions define when a spinor can be written as a swing-twist decomposition 
in respect to a given base vector. 
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